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Abstract 



We study the impact of dimension-six operators on single- and double-Higgs production rates 
via gluon fusion at the Large Hadron CoUider (LHC). If the top- Yukawa couphng is modified 
by some new physics whose scale is of the TeV scale, its effect changes the cross sections of 
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0^' single-Higgs production gg ^ H and double-Higgs production gg HH through the top-loop 

7^ , ■ diagram. In particular, double-Higgs production can receive significant enhancement from the 

, effective top- Yukawa coupling and the new dimension-five coupling ttHH which are induced by 

CO the dimension-six operator. Comparing these results to the forthcoming data at the LHC, one can 

m ' 

■ extract information of the dimension-six operators relevant to the top quark and the Higgs boson. 
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I. INTRODUCTION 



Gauge symmetries of the standard model (SM) have been well confirmed at the CERN 
Large Electron Positron collider (LEP) 111 and the Fermilab TEVATRON. However the 

n n 

mechanism of spontaneous symmetry breaking has not been tested yet [2|, |3| . The vacuum 
expectation value of the Higgs boson triggers electroweak symmetry breaking and generates 
masses of weak gauge bosons, quarks and charged leptons. Search for the Higgs boson is 
the main purpose of the measurement at the CERN Large Hadron Collider (LHC). 

In the SM, coupling constants of the Higgs boson with the weak gauge bosons and the 
matter fields directly relate to their masses. In order to clarify the mass generation mecha- 
nism, an independent determination of the particle masses and their couplings to the Higgs 
boson is important, which will be the subsequent task at the LHC. It will give not only a 
confirmation of the SM but also an indication of new physics beyond the SM. 

At the LHC the dominant production mechanism of the Higgs bosons is gluon fusion 
gg H . The leading contribution to this process comes from top-quark loop diagrams. 
Information of the top- Yukawa coupling can be extracted through this process as a com- 
bination with the Higgs decay branching ratios. The gauge interaction of the Higgs boson 



would be tested through processes of vector boson fusion VV* H{V = W~,Z)^\$ and 

n 

Higgs-strahlung qq' VH [Q]. These processes are promising channels for Higgs searches 
too because of the kinematic advantage in the reconstruction of signals. Measurement of 
the triple-Higgs boson coupling has been discussed in the double-Higgs production mecha- 
nism from gluon fusion, gg HH, at the LHC Q]. In Ref. sl the sensitivity to the Higgs 
boson self-coupling is studied. The authors of this reference conclude that its experimental 
accuracy could reach 20-30% at the SLHC with an integrated luminosity of L = 3000 fb~^ 
for niH = 150-200 GeV. 

Measuring the top- Yukawa coupling accurately is important because the magnitude of 
the coupling constant {yf^ ~ 1) indicates that the physics of top quarks closely would relate 
to that of electroweak symmetry breaking. Lots of models are proposed in this direction j^. 
Measurements of the top- Yukawa coupling would be a key to uncover such possibilities. 
In addition, the Higgs boson self-coupling is of great interest by itself to understand the 
nature of spontaneous symmetry breaking. Its measurement can also be a probe of the new 
physics beyond the SM. The coupling strength is also important being deeply related to the 
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condition of successful electroweak baryogenesis[10|]- 

New physics beyond the SM will be recognized by the discovery of non-SM particles or by 
detecting the deviation from the SM relations between masses and coupling constants. For 
the latter case, the low energy effective theory at the electroweak scale can be described by 
the SM Lagrangian with additional higher dimensional operators. This approach has been 
investigated to analyze non-standard interactions in a model independent way. Leading 
order contributions of such non-standard interactions would be described by the dimension- 



six operators 



111. Constraints on these operators and their phenomenology have been 
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discussed in the literature 

In this paper, we study new physics effects from dimension-six operators on single- and 
double-Higgs production processes, gg H and gg HH . The dimension-six operators 
correct to the top- Yukawa coupling and the triple-Higgs boson coupling, also induce the 
tree level ggH and ggHH vertices. Effects due to the modified top- Yukawa coupling and 
the tree level coupling on the effective ggH vertex are investigated. The former comes from 
color blind new dynamics while the latter can come from some color dependent effects. The 
experimental limits from the LEP precision data and the theoretical bounds such as the 
unitarity bounds are taken into account. We find that the effects on these processes due to 
the dimension-six operators can be significant even under these constraints. In particular the 
double-Higgs production cross section is sensitive to these dimension-six operators. These 
contributions from the dimension-six operators can be distinguished by comparing the data 
for these Higgs boson production channels at the LHC experiments. 

This paper is organized as follows. In Sec. II, we introduce the dimension-six operators 
as a new physics effect. Its experimental and theoretical bounds are discussed. In Sec. 
Ill, numerical evaluations of the effects of dimension-six operators on the Higgs production 
processes at the LHC are shown. Conclusions and discussions are given in Sec. IV. A 
detailed calculation is shown in the appendix. 
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II. EFFECTIVE LAGRANGIAN 



New physics effects on phenomena at the electroweak scale can be described by the higher 



dimensional operators 



The effective Lagrangian is given by 

C, 



cff 



i n>5 



(1) 



where £sm is Lagrangian of the SM, Ci are the coupling strengths of the dimension-n oper- 
ators 0^^\ and A is a cut off scale of the SM. The coefficients of these higher dimensional 
operators can in principle be calculated by assuming new physics models which are defined 
above the scale A. When A is much greater than the electroweak scale, the dimension-six 
operators can give leading contributions to the deviations from the SM. 

If the top- Yukawa interaction is modified by the dimension-six operators, its effect can be 
observed in the processes oi gg ^ H and gg HH . The dimension-six operators relevant 
to the gluon fusion mechanism are 



tG<S> 



G 



jjLUI 



(2) 

(3) 

(4) 



where Q = (u, d)]^, G^^, is the field strength of gluons with SU{3) generators \^{A = 1-8), 
$ is a scalar-iso-doublet with hypercharge F = 1/2, $ = 2x2$*, and v (~ 246 GeV) is the 
vacuum expectation value whose origin may come from color blind dynamics. The effective 
top- Yukawa coupling deviates from the SM value due to On and Oj^^. There is only one 
dimension-six operator Oq 



15| that contributes to ggH and ggHH vertices at the tree level; 



I.e., 



2 / ^i-iu^ ' 



(5) 



whose origin can come from colored new dynamics at the TeV scale. 

The triple-Higgs boson coupling also contributes to double-Higgs production gg HH. 
Dimension-six genuine-Higgs operators change the Higgs self-coupling, which are given 
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by 



C*i = ^5^($^$)5^'($^$), (6) 
0*2 = -^($^$)', (7) 
0^s = {D^^)^^^HD>^^). (8) 

Normalization of the Higgs field is shifted by introduction of the operator O^i. All the Higgs 



interactions are corrected after the wave function renormalization for the Higgs boson [14 1. 
The Higgs potential is modified by C$2; then the triple- Higgs boson coupling is a function of 
C^2 the Higgs boson mass. The coefficient C$2 simply shifts the triple-Higgs coupling. 

The coefficients Ca, C$i and C$2 are free from the current experimental data^. In 
contrast, the coefficient of (9$3 (O^^) contributes to the electroweak rho parameter at the tree 
level (and the one-loop level), which is strongly constrained by the experimental data 13Ul6|. 
The effects of OtG<s? will be measured by top-pair production, gg tt. Here we neglect the 
operators (9^^^, OtG* and (9$3 in the following discussion assuming that these operators 
would be well constrained by the other processes, and we concentrate on the effects due 



to the operators Oti,Oq,i,Oi^2 and Og- In Ref. [15|, the bounds on the operator Oq are 



evaluated through the process gg ^ H ^ WW for > 160 GeV, 



-1.2<a,(^)"<0.5, (9) 



2 

where the scaled couplings are defined as ai = Cij 



A2- 



The theoretical upper bounds on these operators from tree level unitarity 17|] have been 
discussed in the literature Q, [is]]. Such bounds on these less constrained operators are 
given by 

a$i,$2 ^ 47r— . (11) 

If we consider the low energy cut off A = 1-3 TeV, the upper bound on an are 3.0-1.0, 
respectively [l4|. 



^ The coefficient C$i would be determined precisely at the ILC by using the gauge boson association 
processes such as Higgs-strahlung and vector boson fusion. 
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FIG. 1: Feynman diagrams for the single-Higgs production process via gluon fusion. Curled, 
dashed and sohd hnes represent gluons, Higgs bosons and quarks, respectively. Dots denote the 
new physics interaction. In the SM, there is no tree level contact interaction. 
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TABLE I: Parameter sets for the coefficients of the dimension-six operators. 
III. NUMERICAL EVALUATION OF THE HIGGS PRODUCTION PROCESSES 

By introducing the dimension-six operators On and the effective top- Yukawa cou- 
pling is expressed as 

yf = Z^, (^y^ - aa^ , (12) 

where Z$i = (1 + a$i)~^/^. Feynman diagrams for single-Higgs production via gluon fusion 
is depicted in FIG. [H For = 0, the new physics contribution to this process only appears 
in the effective top- Yukawa coupling. Therefore gg ^ H only depends on a^i and a^^ by 
the combination given in Eq. (IT^ . We here define the parameter sets; Set A-Set F, for 
the coefficients of dimension-six operators in TABLE [H Set A corresponds to the SM. The 
effects of Oti are studied by Set B and Set C, and those of On are investigated by Set D 
and Set E. Set F shows those of 0^2- The values of a^i = ±0.5 correspond to the unitarity 
bounds for A = 5 TeV, which change the effective top- Yukawa coupling by about 50%. If 
we take a lower cut off scale (A < 5 TeV), the unitarity constraint on an becomes milder. 

The hadronic production cross section for pp ggX HX is evaluated as a function 
of rrifj at the LHC in FIG. O Detailed calculations are shown in Appendix A. The solid, 
dotted, dashed, long-dashed and dot-dashed curves correspond to the parameter sets Set 
A-Set E, respectively. The peak around mjj ~ 350 GeV in these curves is understood as 
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FIG. 2: The cross section of pp — > ggX HX with ^/s = 14 TeV as a function of the Higgs boson 
mass. Curves denote the cross sections derived in the SM (Set A), and in the SM with anomalous 
dimension-six couplings (Set B-Set E). 

the threshold effect due to the top-quark loop. If the effective top- Yukawa coupling deviates 
from its SM value (Set B and Set C), the cross section can be enhanced by a factor ~ 9/4 or 
suppressed by ~ 1/4 for entire range of the Higgs boson mass. These effects are determined 
only by yf^. The differences from the SM for Set D and Set E are comparable to that for 
top-Higgs coupling for Set B and Set C around ~ 120 GeV. The effects on the cross 
section from Set D and Set E are relatively small compared to Set B and Set C for the larger 
Higgs boson masses. These structures are realized by the interference of the amplitudes 
between the new physics contributions and the SM one. 

In FIGs. |3] and HI we evaluate the statistical sensitivities for anomalous parameters on 
N = La{pp ggX HX)B{H WW,''yy) where the integrated luminosity is assumed 
to be L = 300fb~^. The efficiencies of W bosons and photons are taken as 100% for the 
illustration. Therefore, in these plots, the H —>■ WW decay mode always gives a better 
sensitivity than the ^ 77 decay mode because the decay branching ratios hold the 
relation B{H -> WW) > B{H 77) for > 100 GeV. Note that we do not include 
any backgrounds to calculate the statistical sensitivity. For larger Higgs boson masses, the 
sensitivities become worse due to the decreasing of the cross section. For the WW decay 
mode with rrij^ < 120 GeV, it has also bad sensitivities because of the small branching 
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FIG. 3: The plot of the statistical sensitivity for a^i on N = La{pp ggX — > HX)B{H 
VFVK, 77) where the integrated luminosity is L = 300fb~^. Each curve denotes the la deviation 
from the SM predictions. 




FIG. 4: The plot of the statistical sensitivity for a^. on = La{pp — > ggX HX)B{H 
WW,jj) where the integrated luminosity is L = 300fb~^. Each curve denotes the la deviation 
from the SM predictions. 

fraction of H ^ WW. 

In FIG. m we show the contour plot of the sensitivities in the a^^-aQ plane. There is a 
strong correlation between a^^ and Uq. This means that a^^ can mimic the effect of in 
this process. It is understood by the destructive interference of the top-loop diagram, which 
is shifted by and the tree level diagram, which is induced by the dimension-six operator 
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FIG. 5: The sensitivity plots in the a^-^-aQ plane on = La{pp — > ggX HX)B{H — > l^VF, 77) 
where the integrated luminosity is L = 300fb~^. Each contour represents the la deviation from 
the SM predictions. 






FIG. 6: Feynman diagrams for the double-Higgs production process gg — > HH are depicted. Dots 
represent the new vertices of the dimension-six operators. 

Oq. Therefore, if the deviation from the SM is found in gg H, we cannot distinguish the 
effects of these anomalous couphngs. 

Next, let us discuss the double-Higgs production. In FIG. [HI we show the Feynman 
diagrams for the process gg HH. Invariant amplitudes for the sub-process are given 
in Appendix A. For each diagram, appropriate crossing of external Higgs-boson and gluon 
lines should be taken into account. In the SM, there are only two kinds of topology, the first 
and second diagrams from the left. This process would be used to measure the triple-Higgs 
boson coupling at the SLHC. The vertex function of HHH is modified by the introduction 
of the genuine Higgs operators as 



A 



HHH\^) 



73/2 I y 
^$1 I -^$1 



im^j s + 2m1[ a$2 
2 + 



(13) 



where is the center of mass energy of gg HH. For Set F, the coefficient a^2 is taken 
to be positive to ensure vacuum stability. The value a^2 = +0.5 corresponds to the 140% 
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FIG. 7: The cross section for the Higgs pair production sub-process gg — > HH as a function of 
scattering energy for rrijj = 120 GeV in the SM. The dashed (long-dashed) curve denotes the 
contribution coming from the triangle (box) diagram with the (-|-, -|-) helicity set of gluons. The 
dotted curve is a composition of these two. The thin solid one represents the helicity set (+, — ) in 
the box diagram. 

enhancement of the triple-Higgs boson coupling for = 120 GeV. 

The cross sections for the Higgs boson pair production via the gluon fusion sub-process are 
shown in FIG. [7] in the SM. The dotted and thin solid curves represent the cross section with 
the helicity set (+, +) and (+, — ) of the gluons. For m^j = 120 GeV, the main contribution 
comes from the box diagram (long-dashed) because the triangle diagram (dashed) contains 
the Higgs boson self-coupling, which is proportional to the Higgs boson mass squared, and 
it has the typical behavior of the s-channel process. 

In FIG. [HI we show the effect of the dimension-six top-Higgs interaction on gg HH 
as a function of rrifj for Set B and Set C. These curves are used as the same manner as 
in FIG. El The additional long-dot-dashed curve represents new vertex contribution. The 
SM with gluon helicity sets (+, +) and (+, — ) are also shown in thick and thin solid curves. 
For Set B, it can be seen that the effective top- Yukawa coupling is suppressed in both the 
triangle and the box diagrams. The cross section can also be enhanced by the dimension-five 
interaction vertex of tiHH due to On- On the contrary, for Set C the contributions from the 
box diagrams become large because of the large effective top- Yukawa coupling. The synergy 
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FIG. 8: The cross section for the process gg HH as a function of \/I for = 120 GeV for Set 
B and C. Each curve is given in the same manner as in FIG. [71 The long-dot-dashed curve denotes 
the contribution comes from the new vertex which is induced by the dimension-six operator. The 
SM with hehcity sets (+, +) and (-I-, — ) of the gluons are also shown in thick and thin solid curves 
as a reference. 
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FIG. 9: The cross section of gg HH for rrifj = 120 GeV with a$2 = +0.5 as a function of v s. 
The curves are defined in the same manner as FIG. [HI 

effect of both the contributions enhances the sub-process cross section significantly. 

In FIG. [9l we show the gg HH cross section as a function of \/§ for Set F. Each curve 
is given as in FIG. [H The contribution of '^^^J appears in the triple-Higgs boson vertex. 
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FIG. 10: The cross section of gg — > HH as a function of sub-process energy for rrijj = 120 GeV in 
the SM with Uq = ±0.004. The dotted and long-dashed curves denote the gluon hehcity sets -|-) 
and (-|-, — ) for Set D. The dashed and dot-dashed curves are those for Set E. The SM predictions 
are also given by solid curves. 

Around rrifj ~ 350 GeV, a strong destructive interference between the triangle and the box 
diagrams occurs. The effect in Set F is relatively small as compared to those in Set B and 
Set C. There is also enhancement on the sub-process cross section near the threshold of HH 
production. These effects turn out to give larger contributions to the hadronic cross section. 

For completeness, we also show the case for Set D and Set E in FIG. [TOl The dotted and 
long-dashed curves represent the helicity set +) and (+, — ) of gluons for Set D. Those 
for Set E are given by the dashed and dot-dashed curves. The SM prediction is also shown. 
For Set D, there is a cancellation between the SM contribution and the anomalous tree level 
vertex ggH in the ti threshold region. In Set E, each contribution is constructive in the 
same parameter region. 

Convoluting the CTEQ6M parton distribution function jl9|, the full cross section is 
evaluated for double-Higgs production via gluon fusion. In FIG. [H], we show the hadronic 
cross sections as a function of the Higgs boson mass in the left figure and as a function of the 
dimension-six couplings in the right figure. Significant enhancement from the dimension-six 
operators On and Oq can be seen in both figures for a wide range of parameter space. For 
Set F, the curve and that of the SM one coincide for large m^j values whose structure is 
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FIG. 11: The total cross section of the double-Higgs production pp — > ggX — > HHX as a function 
of the Higgs boson mass in the left figure and the anomalous couplings in the right figure with 
center of mass energy -y/s = 14 TeV. 

easily understood from Eq. f|T3l) . The effect of a$2 is reduced when we take larger Higgs 
boson masses. The contributions from aq,i might be much larger than those of the other 
operators. This effect will be first examined by gauge boson association processes. 

In FIG. [121 we evaluate the statistical sensitivities for the anomalous parameters on 
= La{pp ggX HHX)B{H WW)B{H WW) where the integrated luminosity 
is assumed to be L = 300fb~^. We focus on the anomalous parameters a^l and a^., because 
this process is insensitive to we showed in FIG. [11] Obtained sensitivities are less 

smaller than those in gg ^ H process. However, this process is still sensitive to a^^ for 
~ 150 GeV on some level due to the large enhancement (suppression) in the box 
diagram. On the other hand, it is insensitive to the anomalous parameter a^. 

We show the contour plot of the sensitivities in the a^^-aQ plane in FIG. [131 To compare 
the single- and double-Higgs production processes, we also show the results calculated from 
gg H. By using the insensitivity of gg HH to Qq, the anomalous parameter a^^ can 
be constrained. 

The possibility of measuring the Higgs boson pair production has been discussed in Ref. js] 
to determine the triple Higgs boson coupling constant. Their background analyses can apply 
to our setup because the dominant decay modes of the Higgs boson are almost the same as 
in the SM. In our case, the double-Higgs production process can be observable if the cross 



13 



100 200 300 400 500 600 700 

m„[GeV] 



100 200 300 400 500 600 700 

m„[GeV] 



FIG. 12: The plot of the statistical sensitivity for a^^ on = La{pp ggX HHX)B{H 
WW)B{H — > T^VF) where the integrated luminosity is L = 300fb~^. Each curve denotes the la 
deviation from the SM predictions. 

section receives the above enhancement. However it does not mean the improvement of the 
sensitivity for the triple-Higgs boson coupling. The new vertex tiHH smears the effect of the 
Higgs boson self-coupling. We find that gg HH can still be sensitive to the dimension-six 
top-Higgs interaction. 



IV. CONCLUSIONS AND DISCUSSIONS 



In this paper, we have studied the impact of the dimension-six operators on the Higgs 
production processes via gluon fusion, i.e., gg ^ H and gg HH at the LHC Constraints 
from the current experimental data and the theoretical consistencies on the dimension-six 
operators are taken into account. We find that the contribution from the dimension-six top- 
Higgs operators to single-Higgs production can significantly change the cross section by a 
factor. The double-Higgs production process can also receive a large enhancement from the 
anomalous top-Higgs couplings. The shift of the effective top- Yukawa coupling can enhance 
the cross section significantly. In addition, the new diagrams from the tree level vertex 
tiHH from O^^ result in much larger cross sections than that in the SM. Combined results 
of single- and double-Higgs production can be used to discriminate between the effects of 
dimension-six operators 0^l and Oq. 
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FIG. 13: The sensitivity plots in the a^i-a^ plane on = La{pp — > ggX — > HHX)B{H 
WW)B{H WW) where the integrated luminosity is L = 300fb~^ and = 150 GeV. Each 
contour represents the la deviation from the SM predictions. We also show the contour of the 
sensitivity on = La{pp ggX HX)B{H WW) 

as a reference. 

Finally, we comment on the potential for the coupling measurements at the international 
linear collider. As an optional process the photon-photon collision 77 HH has a similar 
structure to gg HH [20|. This process includes not only the top-quark loop but also 
the W boson loop, so that this kind of enhancement from the top-Higgs interaction may be 
weakened by the W boson loop. The coupling "will also be measured at double-Higgs- 
strahlung e~e^ — > ZHH |2l[|. and W boson fusion e~e^ —>■ vvHH 22|, as well as above the 



photon-photon collision 
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APPENDIX A: CROSS SECTIONS WITH DIMENSION SIX OPERATORS 



We present formulae for the single- and the double-Higgs production processes via the 
gluon fusion mechanism. 

The hadronic cross section at the leading order is calculated by convoluting with the 
parton distribution function (CTEQ6M) as 



6m 



3 ^H^gg' 
H 



T I —g{x)g{T/x), 

X 



(Al) 



where r = m'jj/s, and g{x) is the gluon distribution function in a proton. In the SM with 
higher dimensional operators, the decay rate for the Higgs boson into gluons is calculated 

as 



H-*99 



V2n 



ativ 



V2mtJ m% 



Ami 



H 



Ami 



m 



H 



Co{m 



(A2) 



where Co(s) = Co(0, 0, s, m^, m^, m^) and ^^(/i) is the running strong coupling constant. 



We here use the Passarino-Veltman functions for loop calculations 
corrections are known for the gluon fusion mechanism 



251, 



24| . Large higher order 



26j. The next to the leading 



order (NLO) correction for the cross section can be treated by the K-factor 25|, 



Tv-NLO 
^gg^H 



1 



IT 



11 
Y 



(A3) 



A naive introduction of dimension-six operator breaks renormalizability of the theory. How- 
ever, these operators can be embedded in the more fundamental theory. Thus we here adopt 
the correction only from gluon emissions to the effective ggH vertex in the heavy top-quark 
mass limit as the NLO correction. 

In the effective theory, loop integrations by four momenta are cut off at a some new 
physics scale A, 



(An) 



(27r)4 (P - CY 



-IfC 



A2/C 



t 



1 + ty 



These integrals are related to those in the dimensional regularization (DR) as 



^0 



dx 



dx 



dy 



l-x 



dy 



He (1 + Hc/A^) 



DR 



1 

A^^ 



dz 



1 l + Sffp/A^ 
HI{1 + Hd/A^) 



r)DR_ ^ 
° 2A2^ 



(A4) 

(A5) 
(A6) 
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with 



H' 



c 



( I \2 2 2,2 



E 



D 



(A7) 

[xri + yr2 + zrs) — xr^ — — zr^ + m^, (A8) 

where rj = 'Y^^Pi- The function will appear in the gg — if-ff cross section. We should 
comment on the reduction formulae of the loop integrals in the Passarino-Veltman tech- 
nique. These reduction formulae are fully supported by Lorentz invariance, but the cut off 
regularization generally violates it. In our analysis, we omit the effect of the cut off A in the 
loop integrals. The corrections due to the cut off are small (w/A)^ when the scale A is set to 
be more than 3 TeV, so that the effects of the error turn out to be numerically unimportant. 

Let us discuss the helicity cross section for the sub-process gg HH. The differential 
cross section is calculated as 



daigxgy^HH) \M 



AA 



/ 12 



dt IbTTS^ 

In the SM with dimension six operators (9(i,(9$2 aud Oq, helicity amplitudes M.'^^' are 
given by 



M- 



M — 
Oigtritv 
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(AlO) 
(All) 



where s, t, and u are the Mandelstam variables and the box functions Dx and Dy are 
defined by 

= 4 {2 + 4 Co{s) -m]{s + 2m]j- Sm^) (D^^a ^ ^213 ^ ^132^) 
Am? 



+ 



m 



H 



[2 [i - ml) C* + 2 (m - m|) - {t 



u — m 



H 



) } , (A12) 



Dy = -2 {-s iD^^'^ -su Df^ + 2 s Co(s + 2 (t - m^) C* + 2 (m - m^) 

+ [s-2ml + 8m?) [-2 m? [DI^^ + Df' + Dl'^) + 2 C', 
1 



{-s P D^^^ - s Df^ -s {s-2 mjj) Co{s) 
s {s - 4m|) Q + 2i {i - m]j) Cl + 2u{u- m]j) C^)] ] 



iu — m\j 



(A13) 
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where the scalar loop integrals are listed here: 

(A14) 
(A15) 
(A16) 
(A17) 
(A18) 
(A19) 

The NLO correction for the Higgs boson pair production in the gluon fusion mechanism has 





=Co(0, i, mjj, ml ml mf), 






Co 


=Co(0, u, rriH, ml ml mf), 






Co 


z^/- 2 2 2 2 2\ 

=Co[s, mjj, mfj, m^,m^,m^), 






Dr 


=Do{0, 0, m^, m|, s, t, ml ml 




ml 




=Do(0, 0, m^, mjj, s, u, ml ml 




m^^ 




=Do{0, m|, 0, m^, u, i, ml ml 




2 

m^ 



been studied in the limit of the heavy top-quark mass 27|]. The effect increases the cross 



section to 1.9-2 times that at the leading order. We here take K^^_^^^ = 1.9 as the NLO 
correction to the gg — ^ HH process. 



APPENDIX B: FEYNMAN RULES 

When we introduce the d.„.e,.ion-si. genuine Higg. opetato. 0„ and 0«, the Feyn„.a. 

rules for the Higgs self-interaction are given by [IJ] 

HHH■.-^^v Zl, fz-ll + || E ■ P.. + "-f] (Bl) 

V j<k J 

HHHH : - 6 z 4, [z^l'^ + ^ E ^'^^ " + 2«*2 ) • (B2) 

In the SM with the fermionic dimension-six operator 0^l, the effective top-Higgs interactions 
are obtained as 



ttH:-tZ„[^^-^^^j (B3) 

ttHH:tZl,^, (B4) 

where the wave function renormalization due to the genuine Higgs operator O^i is also taken 
into account. The gluonic operator Oq can induce gluon-Higgs vertices at tree level as 

Gi{Pi)G^{p2)H : - 2z Z^.^p, . p, - ^-^) 6^^ (B5) 

^ V Pi ' P2 / 

G^^{p^)G^ip,)HH : - 2^ Zlfilp, ■ p, - ^^"j 5^^. (B6) 
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